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ABSTRACT 

We present near-infrared (NIR) spectroscopic observations of the blue com- 
pact dwarf (BCD) galaxy Mrk 59, obtained with the TripleSpec spectrograph 
mounted on the 3.5m APO telescope. The NIR spectrum of Mrk 59, which cov- 
ers the 0.90 /im - 2.40 /xm wavelength range, shows atomic hydrogen, molecular 
hydrogen, helium, sulfur and iron emission lines. The NIR data have been sup- 
plemented by a SDSS optical spectrum. We found extinction in the BCD to be 
low [A(y)=0.24 mag] and to be the same in both the optical and NIR ranges. 
The NIR light does not reveal hidden star formation. The H 2 emission comes 
from dense clumps and the H 2 vibrational emission line intensities can be ac- 
counted for by photon excitation. No shock excitation is needed. A CLOUDY 
photoinization model of Mrk 59 reproduces well the observed optical and NIR 
emission line fluxes. There is no need to invoke sources of ionization other than 
stellar radiation. The [Fe ii] 1.257 and 1.643 /im emission lines, often used as 
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supernova shock indicators in low-excitation high-metallicity starburst galaxies, 
cannot play such a role in high-excitation low-metallicity H II regions such as 
Mrk 59. 

Subject headings: galaxies: abundances — galaxies: irregular — galaxies: ISM 
H II regions — galaxies: individual (NGC 4861) — infrared: ISM 



1. INTRODUCTION 



The Blue Compact Dwarf (BCD) gal axy Markarian 59 (Mrk 59) = I Zw 49 belongs 
to the class of cometary BCDs defined by lLoose fc Thuan I (119851 ) as characterized by an 
intense starburst at the end of an elongated low surface brigh t ness stellar body. In the case 
of Mrk 59, the elongated body is the galaxy NGC 4861. krp I Jl966h describes NGC 4861 in 
his Atlas of Peculiar Galaxies as "an object with irregular clumps, resolved into knots with a 
very bright knot (diameter = 1 kpc) at the southeastern end" . The knots are in fact a chain 
of H II regions, most likely resulting from propagating star formation along the galaxy's 
elongated body a nd ending with the h igh-surface-brightness supergiant H H region at the 



southe astern end (INoeske et al. 1 120001 ) . Fabry-Perot maps of NGC 4861 by iThuan et al 



( 119871 ) show that Mrk 59 is ~ 50 times brighter in Ha than any other feature in NGC 4861. 
The Ha velocity field of NGC 4861 shows a northeast-southwest velocity gradient which can 
be interpreted as the rotation of an inclined disk, with the bright end approaching and the 
faint end receding. 



Dinerstein fc Shields I (119861 ) first detected broad Wolf-Rayet stellar feat ures in Mrk 59 



i ndica ting the presence of late nitrogen and early carbon Wolf-R ayet stars. iGuseva et al. 



( 120001 ) . using the spectroscopic observations of 



of Wolf-Rayet stars are present. INoeske et al. 



Izotoy et al. I ( 119971 ) found that several dozens 



(boooh used deep ground-based spectrophoto- 



metric observations of the supergiant H II region to derive an oxygen abundance (log 0/H)j 
= -4.011 ± 0.003 [ or Z m /5, using log (O /H) = -3.3] for the ionized gas, typical of BCDs. 
On the other hand, IThuan et al. I (l2002bl ) have found from FUSE far-UV absorption spectra 
that the oxygen abundance of the neutral gas in Mrk 59 is about a factor of ten lower, (log 
0/H) n = -5.0 ± 0.3. They also did not find any evide nce for the presence o f diffuse molec- 
ular hydrogen in the BCD. IWilcots et al. I (119961 ) and IThuan et al. I (j2002al ) have used the 
VLA to map NGC 4861 and Mrk 59 in the 21 cm line. With a beam size of 15" (780 pc), 
the latter authors found a very clumpy interstellar medium (ISM) with H I column densities 
ranging from several 10 19 to a few 10 21 cm -2 . Despite irregularities, the H I veloc i ty field 
of NGC 4861 resembles that of a rotating disk seen almost edge-on. IThuan et al. I (j2002al ) 



obtained a H I heliocentric velocity of 833 km s 1 for NGC 4861. Then, adopting a Hubble 
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constant of 73 km s" 1 Mpc -1 and correcting for the Virgocentric infall motion, the distance 
of Mrk 59 is 10.7 Mpc (IThuan et al. Il2002al ). Since the apparent magnitude of Mrk 59 is 
rrig = 12.64 (INoeske et al. 1120001 ). then its absolute magnitude M B is equal to -17.51. 



In the near-in f rared (NIR) , NGC 4861 and Mrk 59 have been st udied spectroscopically 
by ICalzetti et al. I (119961) and IVanzi fc Rieke I (119971 ). respectively. ICalzetti et al. I (119961 ) 
(see also ICalzetti Ill997l ) have obtained spectra of NGC 4861 in the wavelength range 1.236 
- 1.340 /mi (J spectrum) and 2.113 - 2.250 /im (K spectrum ). They detecte d the [Fe n] 
1.26 jum, Pa/3 1.28 /im, Br7 2.17 /xm and H 2 2.12 /im lines. IVanzi fe Rieke I (119971 ) have 
obtained a K spectrum of Mrk 59 covering the 2.0 - 2.4 /im wavelength range. Those 
authors detected only three emission lines: the Br7 2.17 /mi hydrogen line, the He I 2.06 
/im line and the H 2 2.12 /tm line, with the last line being barely detected. With advances 
in NIR detector sensitivities, and with the commisionning in 2008 of the new TripleSpec 
spectrograph at the Apache Point Observatory (APO) 3.5-meter telescop^j], which allows 
simultaneous coverage of the JHK NIR bands, we decided to obtain a new NIR spectrum 
of Mrk 59 with considerably more signal-to-noise and wavelength coverage. These improved 
observations will allow us to study in more detail the physical conditions in Mrk 59: its 
extinction in the NIR and the excitation mechanisms of line emission of ionized species in 
the H II region. We can also check for the presence of molecular hydrogen emission. 

We describe the 3.5m APO observations in section 2. In section 3, we discuss the 
extinction and star formation in Mrk 59, both in the optical and NIR ranges, the excitation 
mechanisms for molecular hydrogen and iron emission in the NIR range, and CLOUDY 
modeling of the H II region. We summarize our findings in section 4. 



OBSERVATIONS 



A NIR spectrum of Mrk 59 was obtained with the 3.5 m APO telescope on the night 
of 2008 March 22. The observations were m ade with the TripleS pec spectrograph by the 
instrument commissioning team. TripleSpec (IWilson et al. 1120041 ) is a cross-dispersed NIR 
spectrograph that provides simultaneous continuous wavelength coverage from 0.95 to 2.46 
//m in five spectral orders during a single exposure. A l'. / lx43" slit was used, resulting in a 
resolving power of 3500. 



The A0/AHV standard star HIP 73613 was observed for flux calibration. Spectra of Ar 



1 The Apache Point Observatory 3.5-meter telescope is owned and operated by the Astrophysical Research 
Consortium. 
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comparison arcs were obtained on the same date for wavelength calibration. Since both Mrk 
59 and HIP 73613 are much smaller than the length of the slit, the nod-on-slit technique 
was used to acquire the sky spectrum. Both objects were observed by nodding between 
two positions A and B along the slit. Mrk 59 had a total integration time of 20 minutes. 
During that time the galaxy was observed in four five-minute integrations, alternating in slit 
position following the sequence ABBA. 

The data reduction was performed using IRAF@ and include the following steps. We 
first used the summed frame of the standard star and the IRAF routines IDENTIFY, REI- 
DENTIFY, FITCOORD and TRANSFORM to correct for the distortion of each frame. This 
distortion is an unavoidable consequence of cross- dispersed grating spectrographs. This was 
done for all frames of the standard star, Mrk 59 and the comparison argon lamp. We 
then extracted from each of the distortion-corrected frames the five two-dimensional spec- 
tra corresponding to the five spectral orders. For each order, we use the two-dimensional 
distortion-corrected spectrum of the comparison argon lamp and the IRAF routines IDEN- 
TIFY, REIDENTIFY, FITCOORD, TRANSFORM to perform wavelength calibration and 
correct for the tilt of the slit for each frame. We checked the accuracy of the two-dimensional 
transformation by examining the strong lines in the A and B spectra of Mrk 59. We found 
that the wavelengths of these do not differ by more ~ 0.1 - 0.2A in the A and B spec- 
tra. Then, for each order, the two A and two B distortion-corrected, tilt-corrected and 
wavelength-calibrated frames, each with an exposure time of 5 mn, were summed separately. 
The B+B frame was subtracted from the A+A frame to remove the background. Then, 
the -(B+B) spectrum was shifted along the spatial axis to match the (A+A) spectrum 
and subtracted from it, resulting in the two-dimensional background subtracted spectrum 
(A+A) + (B+B) of Mrk 59 with a total exposure time of 1200 s. The standard star observa- 
tions were processed similarly, resulting in a spectrum with a total exposure time of 300 s. 
Finally, we use the routine CRMEDIAN to remove cosmic ray hits in each spectrum. 

A one-dimensional non-flux-calibrated spectrum of Mrk 59 was extracted from the two- 
dimensional frame using the APALL IRAF routine, within a l'.'lx6" rectangular aperture so 
as to include the brightest star-forming regions. The extraction aperture is shown in Fig. [TJ 
overlaid on an archival HST/WFPC2 image of Mrk 59. We will hereafter refer to it as the 
NIR aperture. 

We simultaneously correct for telluric absorption and flux-calibrate the spectrum of 



2 IRAF is the Image Reduction and Analysis Facility distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities for Research in Astronomy (AURA) under 
cooperative agreement with the National Science Foundation (NSF). 
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Mrk 59 by dividing the non flux-calibrated spectrum of Mrk 59 by that of the standard 
star HIP 73613 and by multiplying the result by the synthetic flux-calibrated spectrum of 
that same standard star. Since there does not exist a published absolute spectral energy 
distribution of the standard star HIP 73613, we have used the synthetic absolute spectral 
energy distribution of the star Vega (a Lyrae) which has a similar spectral type (AOV), and 
scaled it to the brightness of HIP 73613. 

3. RESULTS AND DISCUSSION 

The resulting rest-frame NIR spectra of Mrk 59, corrected for telluric absorption, are 
shown in Fig. [2j Each of the five orders is shown in a separate panel, with the identified 
lines labeled. Each panel shows two spectra. The upper spectrum is the blown-up version, 
corresponding to the flux scale given on the y-axis. It excludes the noisy parts. The lower 
spectrum is the downscaled version of the upper one by a factor of 50. It shows the whole 
wavelength range of the spectrum in each order. 

To have a more complete physical picture of the extinction and the ionization mecha- 
nisms in Mrk 59, we have supplemented the NIR observations with optical ones. We have 
used the optical spectrum of Mrk 59 in the Sloan Digital Sky Survey (SDSS) archives ob- 
tained within a circular aperture of 3" diameter, which covers the wavelength range 0.38 - 
0.93 /im. The optical and NIR apertures are centered on the same spatial location in Mrk 
59 and include the same brightest star- forming regions (see Fig. [1]), so the optical and NIR 
spectra are directly comparable. They overlap in the wavelength range of 0.90 - 0.93 //m, 
allowing us to match one to the other. Specifically, we have adjusted the optical and NIR 
spectra by scaling the latter by a factor 1.03, so that the flux of the [S ill] 0.907 /im emission 
line is the same in both spectra. This procedure works well despite the fact that this line 
is located in a noisy part of the NIR spectrum and its flux is derived with an error of 10% 
(the flux error of the same line in the optical spectrum is smaller, of the order of 3%). The 
combined redshift-corrected optical and NIR spectrum is shown in Fig. [3J Gaps in the NIR 
spectrum are regions of strong absorption by the telluric lines. 

Emission-line fluxes in both the optical and NIR ranges were measured by using the 
SPLOT routine in IRAF. The errors of the line fluxes are based on the photon statistics in 
the non-flux-calibrated spectra. Additionally, we have also included an error of 10% to all 
lines. This error is that of the [S ill] 0.907 /xm emission line flux, and including it takes into 
account the fact that this line was used to adjust the NIR spectrum to the optical one. The 
observed fluxes F(X) of emission lines derived from the adjusted optical and NIR spectra 
are shown in the second columns of Tables [T] and [2J They are normalized to the observed 
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H/3 flux, F(R/3) = (2.305±0.008)xl(T 13 erg s" 1 cm- 2 , and multiplied by a factor of 100. 
In Table we have indicated the emission lines that are affected by telluric absorption. 
Note that several important emission lines in the optical SDSS spectrum are absent: [O n] 
0.373 /mi is out of the wavelength range of the SDSS spectrum, while the Ha 0.656 /mi and 
[O ill] 0.501 /mi emission lines are clipped. For these lines, we have adopted the observed 
fluxes (relative to H/3) derived by llzotov et al. I (119971 ). This is a vali d procedure because 



the relative intensities of the common emission lines in the SDSS and llzotov et al. 
spectra are in good agreement within the errors. 



(1997) 



Previously, Mrk 59 has been observed in the NIR by IVanzi &: Rieke I (119971 ). The line 
fluxes of the three emission lines they detected (Br7, He I 2.06 /im an d H2 2.12 /mi) ar e lower 
than ours by a factor of 2.5 - 3.0. However, the continuum level in the lVanzi &: Rieke I (119971 ) 
spectrum is a factor of two higher than that in our spe ctrum. While th e high er continuum 
level can be explained by the larger aperture used by IVanzi fe Rieke I (119971 ) (2'.'4xl2" as 
compared to our l'/lx6" extraction aperture), aperture effe cts cannot account for the lower 



line fluxes. Th e reasons for such differences are not clear. ICalzetti et al. I (119961 ) (see also 



Calzetti I Il997l ) have also observed the Br7 and [Fe 11] 1.257 /mi lines in Mrk 59. Their 
considerably larger aperture (6 / .'9x20 // ) gives continuum and emission line fluxes that are 4-5 
times larger than our fluxes. However, the Br7 to [Fe 11] 1.257 /xm line r atio is consistent with 
ours within the errors. We obtain F(Br7) /F([Fe ll]1.257) = 10±5 while ICalzetti I (119971 ) gets 
15±6. 

The presence of many emission lines in the optical and NIR spectra (Figs. [U [3] and 
Tables [1] and [2]) allows us to study extinction and star formation, H2 and Fe II emission and 
excitation mechanisms in Mrk 59. We discuss these issues in the following sections. 



3.1. Extinction and star formation 

The fact that the NIR spectrum of Mrk 59 has been obtained simultaneously over the 
entire JHK wavelength range and that it is possible to directly match the optical spectrum 
to it, using common emission line fluxes in both spectra, eliminates uncertainties introduced 
by the use of different apertures in the optical and NIR observations and by the adjusting of 
the continuum levels of NIR spectra obtained separately in different orders. In all previous 
studie s of BCDs, the N IR spectra of BCDs were obtained in separate JHK observations 
(e.g., IVanzi et al. 1120021 ) . and there was no wavelength overlap between the optical and NIR 
spectra. The elimination of these adjusting uncertainties permits us to compare directly the 
optical and NIR extinctions. 
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In the second column of Table [3] we show the observed fluxes F(X) of the strongest 
hydrogen lines in both optical and NIR spectra which are not affected by telluric absorption. 
For comparis on, the fourth column of T able [3] shows theoretical recombination flux ratios, 
calculated by iHummer Storey I (119871 ) for case B with an electron temperature T e = 15000 



K and an electron numbe r density N K = 10 cm 3 , which approximate well the values 



observed in Mrk 59 (e.g., Ilzotov et al. 1119971 ). It is seen from Table [2] that the observed 



and theoretical ratios ar e close to each other for all hydrogen lines, suggesting very little 
extinction. Adopting the lCardelli et al. I (119891 ) reddening curve with R(V) = A(V)/E(B — 
V) = 3.2, we derive an extinction C(H/3) = 0.11 from the optical spectrum (Tabled]). This 
value of C(H/3) c orresponds to MV) = 0.24 mag, consistent with the value A(V) = 0.2 



mag obtained by Ilzotov et al. I (119971 ). The extinct ion- corrected fluxes /(A) of the strong 
hydrogen lines are shown in the third column of Table [3J There is a very good agreement 
between the corrected and theoretical recombination values of the line fluxes in both the 
optical and NIR ranges. This implies that a single C(H/3) [or -A(V)] can be used over the 
whole 0.38 - 2.40 /mi range to correct line fluxes for extinction. The extinction-corrected 
fluxes J (A) of all emission lines are shown in the third column of Tables [1] and [2] respectively 
for the optical and NIR ranges. The fact that A(V) does not increase when going from the 
optical to the NIR wavelength ranges implies that the NIR light does not probe more extinct 
regions with hidden star for mation as compared to t he optical light. This appears to be a 



general result for BCDs (e.g. lVanzi et al. Il2000u2002l ). Even in the extreme case of the BCD 
SBS 0335-052E where far-infrared observations imply a very large extinction — 15-2 



magi and three times as much hidden star formation as seen in the visible ( IThuan et al 



19991 ; iHouck et al. 1120041 ). the A(V ) = 0.73 mag deriv ed in the NIR is only slightly greater 
than the optical A(V) = 0.55 mag JVanzi et al. Iboooh . 



Using the equivalent width of 154 A of the H/3 line, Starburst99 ( Leitherer et al. Ill999l ) 
gives a cluster age of a bout 4.5 Myr in th e instantaneous burst case, in good agreement 
with the a ge derived by [r3arth et al. ( 1994 ) for Mrk 59. The presence of Wolf-Rayet stars 



in Mrk 59 ( Guseva et al. 200ol ) also implies that the age of the youngest stellar clusters 



is 

greater than 3 M yr. Another age indica tor is the clear presence of the COA2.3 / um absorption 
feature. Models (ILeitherer et al. Ill999l ) predict the appearance of this feature only after 10 
Myr. This implies the presence of slightly older clusters in Mrk 59 which do not contribute 
substantially to the ionization of the gas. 

The total Ha luminosity of Mrk 59 correcte d for both Galactic and internal extinction 
at a distance of 10.7 Mpc is 7.4 xlO ergs s" - 1 rtGil de Paz et al. Il2003l l This co rresponds 
to a s tar formation rate (SFR) of stars more massive than 5 M & of 0.017 M & yr _1 (ICondon 



19921 ). The total SFR including all stars more massive than 0.1 M Q is 0.09 M yr _1 , typical 
of SFRs of BCDs. The thermal radio continuum flux at 1.4 GHz corresponding to the 
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extinction-corrected H/3 flux is expected to be about 6 mJy (jCondon II1992T ). iThuan et al. 



( j2002al ) have observed a 1.4 GHz flux of 10 mJy in Mrk 59. Thus most of the radio continuum 
emission is of thermal nature. The remaining flux probably comes from supernova remnants 
or some hidden star formation. We fa vor the latter hypothesis as the spectral index is small, 
equal to —0.15 (IThuan et al. Il2002al ). while that of supernova remnants is larger, ~ —0.8. 



3.2. H 2 emission 



Molecular hydrogen lines do not originate in the H II region, but in neutral molecular 
clouds. In the near-infrared, they are excited through two mechanisms. The first one is the 
thermal mechanism consisting of collisions between neutral species (e.g., H, H 2 ) resulting 
from large-scale shock s driven by powerful stell ar winds, supernova remnants or molecular 



cloud collisions (e.g. IShull fc Beckwith Ill982l ). The seco nd one is the fluorescent mech- 



anism due to the absorption of ultraviolet photons (e.g. iBlack fc Dalgarno Ill976l ). It is 
known that these two mechanisms mostly excite different roto- vibrational levels of H 2 . By 
co mparing the observed line ratio s with those predicted by models, such as those calculated 
by IBlack &: van Dishoeck I (119871 ). it is possible to discriminate between the two processes. 
In particular, line emission from the vibrational level v—2 is virtually absent in collisionally 
excited spectra, while they are relatively strong in fluorescent spectra. 

We detect four H 2 emission lines in our NIR spectrum: the 2.034 /j,m 1-0 S(2), the 
2.122 /an 1-0 S(l), the 2.223 fjm 1-0 S(0) and the 2.248 /mi 2-1 S(l) lines. The 2.034/2.122, 
2.223/2.122 and 2.248/2.122 line flux ratios are close to 0.5 (Table EJ, in agreement with the 
values expected for fluorescent excitation, but si gnificantly higher than those predicted for 
thermal excitation (IBlack fc van Dishoeck 1119871 ). In particular, the presence of the 2.248 
/tm 2-1 S(l) line strongly favors the fluoresc ent mechanism. We conclude that the H 2 lines in 
Mrk 59 are mainly excited by fluorescence. IVanzi et al. I (120081 ) also found the fluorescence 
excitation process to be responsible for th e presence of H ? lin es in II Zw 40. This galaxy 
has an equivalent width EW(H/3) = 258A ( Guseva et al. II2006I ). implying that its starburst 
has a slightly younger age of ~ 3 Myr. How do w e reconc ile the fact that H 2 emission lines 
are seen in the NIR spectrum while IThuan et al. I (l2002bl ) did not detect any H 2 absorption 
line in their FU SE spectrum of Mrk 59? The situation is similar to the one in the BCD 
SBS 0335-052E (IThuan et al. 1120051 ) where H 2 emission is detected in the NIR, but no UV 
absorption lines are seen. The explanation for the lack of UV absorption lines in Mrk 59 
is the same as for SBS 0335-052E: the NIR H 2 emission does not come from a uniform 
low-density neutral interstellar medium but from dense clumps. FUSE observations are not 
sensitive to such a clumpy H 2 distribution. They can only probe diffuse H 2 along the lines of 
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sight to the several hundreds of UV-bright massive stars in Mrk 59. As in the case of other 
BCDs, the absence of diffuse H 2 in Mrk 59 can be explained by the combined effects of a 
low H I density, a scarcity of dust grains on which H 2 molecules can form, and a large UV 



flux that destroys molecules (IThuan et al. Il2002bl ). 



3.3. Stellar photoionization CLOUDY model of the H II region 



We next examine the excitation mechanisms for the NIR emission lines in the H II 
region of Mrk 59. For this purpose, we have c onstructed a st ellar photoionization mode l 



for the H II region, using the CLOUDY code of iFerland I (119961 ) and iFerland et al. I (119981 ) . 



We list in Table H] the input parameters of this model. Q(R) is the log of the number of 
i onizing photons per second, calculated from the FL3 luminosity at a distance of 10.7 Mpc 
( IThuan et al. Il2002al ). T e ff is the effective temperature of stellar radiation, N e is the electron 
number density and / is the filling factor. The remaining input parameters are the number 
ratios of different species to hydrogen. We varied the abundances of different elements to 
obtain agreement between the observed and predicted line fluxes. Exceptions are the oxygen 
abundance for which we use the value derived from the SDSS spectrum, and the carbon and 
silicon abundances. There are no strong emission lines o f these species in the opti c al and NIR 
ranges. We have therefore taken their abundances from iKobulnickv fc Skillman ( 1998 ). For 
the ionizing stellar radiation, we have adopted a iKurucz I Jl979T stellar atmosphere model 
with T e //=50000 K. We have also experimented wit h other forms of the i onizing stellar 
radiation such as those given by Starburst99 models (ILeitherer et al. 1119991 ) or by models 
based on Geneva ste llar evolutionary tracks with a heavy element mass fraction Z = 0.004 
(IMeynet et al. II 19941 ). In all cases, we were able to reproduce well the observed line intensities 
by varying only slightly the input parameters. Thus our results are robust and do not depend 
sensitively on a particular adopted grid of models. 

The CLOUDY-predicted fluxes of the optical emission lines are shown in the fourth 
column of Tabled! Comparison of the third and fourth columns of that table shows that, in 
general, the agreement between the optical observations and the CLOUDY predictions is very 
good. The most notable difference is for sulfur. The model reproduces well the [S ill] emission 
line fluxes, but the predicted fluxes for the [S n] nebular lines are lower by a factor of ~ 2 than 
the observed ones. A probable cause for this discrepancy is the poorly known dielectronic 
recombination coef ficients of sulfur, on which the ionic abundances depend sensitively (e.g., 



Izotov et al. 1120061 ). 



Examination of the third and fourth columns of Table [2] shows that the same CLOUDY 
model reproduces also very well the observed emission-line fluxes in the NIR spectrum. This 
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agreement implies that a H II region model including only stellar photoionization as ionizing 
source is able to account for the observed fluxes both in the optical and NIR ranges. No 
additional excitation mechanisms such as shocks from stellar winds and supern ova remnants 



are ne eded. This situation is again similar to that in the BCD II Zw 40 where IVanzi et al. 
( 120081 ) found that the interstellar medium (ISM) is mainly photoexcited by stars. 



3.4. [Fe n] line emission 



The [Fe n] 1.257 and 1.643 fim emission lines are clearly detected (Fig. [T]). They have 
often been used as indicato r s of the importance of sho ck excitation relative to photoionization 
( jMoorwood fe Oliva Ill988l ; lOliva et al. Ill990l . 1200 ll ). In the Galaxy, the intensities of these 
lines are considerably smaller in regions where photoionization dominates than in those where 
shock excitation dominates, such as in SN remnants. For example, the ratio of the [Fe ll] 1.643 
i ttm line to Br7 is 0.06 in the Orion nebula, w hile it is more than 30 in Galactic SN remnants 
( jMoorwood fe Oliva Ill988u01iva et al. Ill990l ). This is because fast shocks propagating in the 
diffuse ISM can destroy dust grains through sputtering processes or grain-grain collisions, and 
replenish the ISM with gas-phase iron. The gaseous iron is then collisionally excited in the 
cooling post-shock gas, and produces the observed infrared emission. 

In the Mrk 59 spectrum, the ratios of the [Fe ll] 1.257 fim and [Fe ll] 1.64 3 jum line 



i ntens ities to that of Br7 have the very low values of 0.099 and 0.118, respectively. ICalzetti 



( 119971 ) also found that the [Fe ll] 1.257 /im line flux is some 30 times smaller in Mrk 59 (for a 
given Pa/3 flux) as compared to the remaining galaxies in her sample of 19 starburst galaxies, 
all having higher metallicity (about solar) and lower excitation. Does this mean that there 
is no significant SN shock excitation in Mrk 59? If this is the case, this raises a puzzle since 
the age of the ionizing cluster in Mrk 59 is 4.5 Myr. For a Salpeter IMF and an upper mass 



limit of 100 M , Starburst99 (ILeitherer et al. 111999 ) gives a SN rate of 1.4xl0~ 4 yr^ 1 . Since 
the first SN explodes around 3 Myr, we would expect ~ 200 SNe to have occurred by the 
age of 4.5 Myr. 

Examination of Table [2] shows that the stellar photoionization CLOUDY model repro- 
duces reasonably well the observed [Fe ll] 1.257 and 1.643 ^m line fluxes. We have also run 
a CLOUDY model where the ionization radiation does not come from stars like in the model 



discussed in Section 3.3, but from fast radiative shocks, as modeled by lAllen et al. I (120081 ). 



We use the slowest available shock model from I Allen et al. I (120081 ). with a velocity of 100 
km s _1 . With the observed log Q = 51.953 (Table SJ), we again obtain [Fe ll] 1.257 //m/Br7 
and 1.643 //m/Br7 line flux ratios of ~ 0.1, similar to those obtained with the model with 
stellar ionizing radiation (Table |2J). It is clear that, in a high-excitation H II region as Mrk 
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59, the [Fe n] 1.257 and 1.643 /im emission lines are not particularly sensitive to the presence 
of shocks, since CLOUDY models with a stellar ionizing spectrum or with a shock ionizing 
spectrum produce the same line intensities. However, the shock model fails to reproduce 
the strong forbidden emission line fluxes, for example overpredicting the [O n] 0.373 /im 
emission line by 40% and underpredicting the [O ill] 0.501 /im emisssion line by 80%. It is 
clear that shocks, although likely present in Mrk 59, cannot be the main source of ionization 
in it. Stellar radiation is. 

We thus conclude that the [Fe n] 1.257 and 1.643 jttm emission lines cannot be used as 
shock indicators in hi gh-excitation low-meta llicity H II regions. This means, for example, 



that the Starburst99 (ILeitherer et al. Ill999l ) predictions for the [Fe n] 1.257 line flux as a 



function of the supernova rate cannot be ap plied to M r k 59. This is because these predictions 
are based on mean relations established by Calzettil Jl997l) fo r [Fe n] line emission in low- 
excitation high-metallicity starburst galaxies. ICalzetti I (119971 ) found that the [Fe n] 1.257 
line flux of the high-excitation low-metallicity BCD Mrk 59 (= NGC 4861) deviates strongly 
from that mean relation, being a factor of 30 below the [Fe n] 1.257 line flux - Pa/5 line 
flux mean relation. This could be due to two reasons: 1) the lower metallicity in Mrk 59 
results in a lower Fe abundance; and 2) the higher excitation of Mrk 59 results in an iron 
ionic composition which is very different from that for other galaxies. Our CLOUDY stellar 
photoionization model predicts that most of the iron in Mrk 59 is in the Fe IV ionization 
stage, there being about 10 times more Fe IV ions than Fe II ions (by volume), and about 3 
times more Fe ill ions than Fe II ions. 



4. CONCLUSIONS 

We discuss here the results of near-infrared (NIR) spectroscopic observations of the 
blue compact dwarf (BCD) galaxy Mrk 59. The NIR data have been supplemented by the 
optical spectrum from the Sloan Digital Sky Survey (SDSS), resulting in a total wavelength 
coverage of 0.38 - 2.46 /im. The overlap in the 0.90 - 0.93 /im wavelength range allows to 
adjust the flux level of the NIR spectrum to that of the optical spectrum, using fluxes of 
common emission lines. 

We have arrived at the following conclusions: 

1. Using the Balmer decrement of hydrogen emission lines in the optical range, we 
find that the extinction in Mrk 59 is low, amounting to A(V) = 0.24 mag. Correcting the 
observed hydrogen emission line fluxes in the NIR range for this optically-derived reddening 
gives values that are very close to the theoretical recombination values. Thus, extinction in 
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Mrk 59 is the same in both optical and NIR ranges. The NIR light does not probe more 
extinct regions with hidden star formation as compared to the optical light. 

2. We have detected four molecular hydrogen emission lines in the NIR spectrum of Mrk 
59. Comparison of the observed fluxes with modeled ones suggests that the main excitation 
mechanism of H2 emission in Mrk 59 is fluorescence. The H2 emission comes not from a 
uniform low-density neutral medium but from dense clumps. 

3. A CLOUDY model with a pure stellar ionizing radiation reproduces well the observed 
emission line fluxes in both optical and NIR ranges. Shocks are likely present in the H II 
region, but play a minor role in the ionization. 

4. CLOUDY stellar ionizing and shock excitation models show that the [Fe 11] 1.257 and 
1.643 /iin emission lines, often used as SN shock indicators in low-excitation high-metallicity 
starburst galaxies, cannot play such a role in high-excitation low-metallicity H II regions 
such as Mrk 59. Both types of models predict similarly low [Fe 11] emission line fluxes for 
the latter category of objects, so that these line fluxes cannot discriminate between the two 
mechanisms. 

We thank The TripleSpec team for their superb work in building TripleSpec and the 
APO staff for their valuable assistance during commissioning. Y.I.I, is grateful to the staff of 
the Astronomy Department at the University of Virginia for their warm hospitality. Support 
for this work is provided by NASA through contract 1263707 issued by JPL/Caltech. 
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Table I. Fluxes of optical emission lines 
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Table 1 — Continued 



Ion 


100xF(A)/F(H/3) a 


100x/(A)//(H/?) b CLOUDY 


0.655 [N II] 


1.70±0.09 


1.58±0.08 


1.58 


0.656 Ra c 


307.20±8.77 


285.15±8.71 


283.73 


0.658 [N II] 


4.68±0.16 


4.34±0.16 


4.65 


0.668 He I 


3.06±0.12 


2.83±0.12 


2.74 


0.672 [S II] 


9.85±0.30 


9.09±0.30 


4.74 


0.673 [S II] 


7.22±0.23 


6.67±0.23 


3.59 


0.707 He I 


2.52±0.10 


2.30±0.10 


3.15 


0.714 [Ar ill] 
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0.29 


0.850 Pal6 
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0.34 


0.855 Pal5 


0.47±0.15 


0.41±0.14 


0.45 


0.860 Pal4 
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0.99 


0.867 Pal3 


0.73±0.15 


0.63±0.14 


1.03 


0.875 Pal2 
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1.20 


0.886 Pall 


1.46±0.17 


1.26±0.16 


1.49 


0.901 PalO 


2.08±0.19 


1.79±0.18 


1.92 


0.907 [S ill] 


16.96±0.50 


14.62±0.54 


15.06 



a F(H/?) = (2.305±0.008)xlO- 13 erg s" 1 cm" 2 . 
b Corrected for extinction with C(H/5)=0.11. 
c Observed flux is from llzotov et al. I (119971 ). 



Table 2. Fluxes of near-infrared emission lines 



Ion 100xF(A)/F(H/?) a lOOx J(A)/J(H/3) b CLOUDY 



0.907 [S m] 


16.98±1.63 


14.63±1.40 


15.06 


0.953 [S m] c 


43.43±4.27 


37.13±3.88 


37.35 
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3.26±0.43 


2.79±0.41 


3.66 


1.005 Pad 
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5.47 
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0.15 
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1.257 [Fe II] 
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1.279 He I 
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0.17 
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0.22 


1.588 Brl4 


0.35±0.07 
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0.47 


1.611 Brl3 
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0.49 


1.641 Brl2 
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0.43±0.08 


0.58 


1.643 [Fe n] 


0.36±0.07 


0.30±0.07 


0.33 


1.681 Brll 


0.78±0.11 


0.63±0.10 


0.71 


1.700 He I 


0.38±0.09 


0.30±0.07 


0.28 


1.732 He I 


0.10±0.04 


0.08±0.03 




1.736 BrlO 


1.02±0.13 


0.82±0.09 


0.92 


1.745 [Fe ll]+He I 


0.07±0.05 


0.06±0.04 




1.876 Paa c 


33.63±3.18 


26.97±2.78 


31.74 


1.944 Bt5 c 


2.16±0.37 


1.73±0.28 


1.77 


2.034 H 2 1-0 S(2) 


0.13±0.04 


0.10±0.04 




2.058 He l c 


1.32±0.16 


1.07±0.14 


1.45 
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Table 2 — Continued 



Ion 


100xF(A)/F(H/3) a 


100x/(A)//(H/9) b CLOUDY 


2.122 H 2 1-0 S(l) 


0.26±0.09 


0.21±0.08 


2.165 Br 7 


3.18±0.32 


2.54±0.30 2.64 


2.223 H 2 1-0 S(0) 


0.10±0.06 


0.08±0.05 


2.248 H 2 2-1 S(l) 


0.15±0.06 


0.11±0.05 



a F(H/3) = (2.305±0.008)xl0- 13 erg s" 1 cm- 2 . 
b Corrected for extinction with C(H/3)=0.11. 
c The observed flux is affected by telluric absorption. 
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Table 3. Fluxes of strong hydrogen lines 



Line 100xF(A)/F(H/3) a 100x/(A)//(H/3) b Theor. c 





a) Optical spectrum 






0.410 m 


24.7 


26.2 


26.2 


0.434 H 7 


46.8 


48.7 


47.3 


0.434 H/3 


100.0 


100.0 


100.0 


0.875 Pal2 


1.1 


1.0 


1.0 


0.886 Pall 


1.5 


1.3 


1.3 


0.901 PalO 


2.1 


1.8 


1.8 




b) NIR spectrum 






1.094 Pa7 


9.1 


7.6 


8.6 


1.282 Pa/3 


16.4 


14.6 


15.2 


2.165 Br 7 


3.2 


2.5 


2.5 



a F(H/3)=2.305xlO" 13 ergs s^cm" 2 . 
b Corrected for extinction with C(H/3)=0.11. 
c Recombination ratios for 15000 K and N e = 100 cm -3 



from iHummer &: Storey I (Il987l ) 
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Table 4. Input parameters for stellar photoionization CLOUDY model 



Parameter 


Value 


Q(H) 


51.953 


Teff, K 


50000 


A^ e , cm -3 


100 


/ 


0.07 


log He/H 


-1.07 


log C/H 


-4.48 


log N/H 


-5.34 


log O/H 


-3.97 


log Ne/H 


-4.72 


log Si/H 


-5.35 


log S/H 


-5.69 


log Cl/H 


-7.59 


log Ar/H 


-6.32 


log Fe/H 


-6.22 
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Mrk 59 
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Fig. 1. — The extraction aperture of l'.'lx6" for the one-dimensional NIR spectrum is shown 
by a rectangle. The 3" diameter circular aperture of the SDSS optical spectrum is also 
shown. Both NIR and optical spectra are overlayed on an archival HST/WFPC2 image of 
Mrk 59. North is up and East is to the left. 
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Fig. 2. — The 3.5m APO/TripleSpec NIR spectrum of Mrk 59 in five orders. In each 
panel, the noisy regions of the upper spectrum are omitted. They are caused by insufficient 
sensitivity or strong telluric absorption. The flux scale on the y-axis corresponds to the 
upper spectrum. The lower spectrum is downscaled by a factor of 50 as compared to the 
upper spectrum. It is shown for the whole wavelength range in each order. 
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Mrk 59 
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Fig. 3. — Spectral energy distribution of Mrk 59 in the wavelength range 0.39 yum - 2.4 
yum. The optical spectrum in the wavelength range 0.39 /xm - 0.93 yum is from the SDSS. 
The other spectra have been obtained with the 3.5m APO telescope and the TripleSpec 
spectrograph. 



